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Abstract 

A model of two consumer-resource systems linked by interspecific interference 
competition of consumers is considered. The basic assumption of the model is 
that the dynamics of the resource is much slower than that of the consumer. 
In the absence of interaction each consumer-resource pair has a unique sta¬ 
ble steady state which is completely nonoscillatory. When weakly coupled, 
the consumer-resource pairs are shown to exhibit sustained low-frequency syn¬ 
chronous antiphase relaxation oscillations. 
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1. Introduction 


Recently, synchronization of coupled oscillators has been a subject of exten¬ 
sive study, not only due to the ubiquity of this phenomenon, but also owing to 
importance of its applications in engineering and biology [l|, 12;, 3; \% |5|]. The con¬ 
ventional assumption of the theory of synchronization is that in the uncoupled 
state each elementary unit of the linked system is oscillatory. However no less 
interesting are the systems where coupling is essential for the very emergence 
of oscillations and not only for their synchronization and phase adjustment. 
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The first example of coupling-induced periodicity has been propounded by 
Smale [f|. His abstract model of a biological cell involves chemical kinetics of 
four metabolites, such that the reaction equations for the set of metabolites have 
a globally stable equilibrium. The cell is “dead”, in that the concentrations of 
its metabolites always relax to the same fixed levels. When two such cells are 
coupled by linear diffusion terms, however, the resulting equations are shown to 
have a globally stable limit cycle. The concentrations of the metabolites begin 
to oscillate, and the system becomes “alive”. 

Since that time, triggered by Smale’s seminal work, a number of plausible 
models have been proposed in which coupling of identical nonoscillating cells of 
concrete nature could generate synchronous oscillations. The majority of these 


models concern neural cells with excitable membrane 


Chua 


HQ, fl. 


Szatmari and 


suggested an apt term “awakening dynamics” for the phenomenon. 

The subject of the present paper is an emergence of collective oscillations 
in a simple system of two coupled nonoscillatory consumer-resource pairs. Our 
choice of coupled consumer-resource equations as a matter of enquiry is dic¬ 
tated primarily by the abundance and importance of consumer-resource rela¬ 
tions. Consumer-resource communities are the building bricks of ecosystems. 
Depending on a specific nature of the involved consumer-resource interactions, 
they can take the forms of predator-prey, herbivore-plant, parasite-host, and 
exploiter-victim systems flo| . However applications of the consumer-resource 
models extend far beyond the ecology and are found wherever one can speak 
of win-loss interactions. In its broad meaning, resource is any substance which 
can lead to increased growth rate of the consumer as its availability in the en¬ 
vironment is increased. As this takes place, the resource is certainly consumed. 
Consuming the resource means tending to reduce its availability. When care¬ 
fully examined, consumer-resource models are identified in the following fields: 


epidemiology 
and electrons 

H,p. 


H ected ar 
ch. 6]), 


and susceptible [ll|, ch. 10]), laser dynamics (photons 


labor economics (share of labor and employment rate 
28]), theoretical immunology (antigens and B lymphocytes 114 p. 299]), 
kinetics of chain chemical reactions (free radicals and lipid molecules jl5|), and 
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in numerous other studies from diverse disciplines. 


So far as we know, examples of coupling-induced synchronization of intrin¬ 
sically nonoscillatory populations have never been proposed. We are going to 
show that interaction in a form of density-dependent cross-losses may drive two 
nonoscillatory consumer-resource pairs into synchronous periodic pulsing. 

2. The model 

Of all types of interactions between individuals of the same population (in¬ 
traspecific interactions) or individuals of different populations (interspecific in¬ 
teractions) of the same trophic level competition is most commonly encountered. 
In a broad sense, competition takes place when each species (individual) has an 
inhibiting effect on the growth of the other species (individual). An inhibiting 
effect should be understood to mean either an increase in the death rate or a 
decrease in the birth rate. 

Consider the famous consumer-resource equations proposed by MacArthur 




(la) 



(lb) 


i = 1,..., n. 


Here dots indicate differentiation with respect to time t, x t represents the total 
biomass of jth resource (prey), yi stands for the total biomass of ?’th consumer 
(predator) species, the constant rj defines the growth rate of jth resource, Kj 
is the carrying capacity of jth resource, Cy is the rate of uptake of a unit of jth 
resource by each individual of ith consumer population, wj 1 is the conversion 
efficiency parameter representing an amount of jth resource an individual of ith 
consumer population must consume in order to produce a single new individual 
of that species, bi is the loss rate of *th consumer due to either natural death 
or emigration. All parameters in Jl]) are nonnegative. 

MacArthur H assumed population dynamics of the resources to be much 
faster than that of the consumers which enabled him to approximate Xj in (flbl) 
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by its quasi-steady-state value derived by setting the right-hand side of Hal) to 
zero. As a result, he succeeded in reducing slow-scale subsystem of equations 
(flbl) to the well-known Lotka-Volterra-Gause (LVG) model fill 



where a is = c ii c sj( w j K j! r j) anci = Djli CijWjKj-bi (i,s = 1,..., n). 
Resources do not enter LVG equations explicitly being parameterized by carry¬ 
ing capacities. 

More recently, such an asymptotic reduction has also been carried out for a 
model of competition where species (with continuous trait) consume the com¬ 
mon resource that is constantly supplied, under the assumption of a very fast 
dynamics for the supply of the resource and a fast dynamics for death and 
uptake rates [h| . 

Consumer-resource model © assumes that competition between consumer 
species is purely exploitative: individuals and populations interact through uti¬ 
lizing (or occupying) a common resource that is in short supply. Quite on the 
contrary, LVG model © describes competition strictly phenomenologically, as 
direct interference where consumers experience harm attributed to their mutual 
presence in a habitat (e.g. through aggressive behavior). However we should 
stress that neither does MacArthur’s reduction claim that interference compe¬ 
tition entirely results from “more fundamental” trophic competition, nor does 
it urge us to hastily consider direct competition as some “derived” concept. 
What it states is that when the dynamics of the consumers are associated with 
a slow time scale, the effects of exploitation competition are indistinguishable 
from those of interference competition. And at slow-time scale, coefficients a,i S 
of © merely would add to “true” interference coefficients a' is if the interference 
is accounted for properly in (US- 

Most mathematical models dealing with coupled consumer-resource pairs or 
multilevel trophic chains ignore contributions of intraspecific and interspecific 
interference effects in consumers. Indeed, the empirical data like Q do indicate 
that a ’ i3 may be negligible in comparison with ai s . Nevertheless, literature ad- 
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vocating the explicit accounting for direct interference shows that incorporation 
of self-limitation and cross-limitation terms in the equations at the consumere’ 
level can provide for the stable coexistence of many species on few resources 


M P- 31; 


21 


231. 


Moreover, if we are to assume dynamics of the resources to be much slower 
than that of the consumers, it is likely that we have to introduce interference 
competition terms in subsystem m- 

Consider the following modification of © representing coupled two-consumer, 
two-resource equations: 


Xi =P1 - {ciyi + q\)xi, 

(3a) 

X2 = P2 - (c 2 J /2 + q2)x 2 , 

(3b) 

yi = (cxWiXi - bi - diyx - h 2 y 2 )yi, 

(3c) 

2/2 = ( c 2 w 2 x 2 -b 2 - d 2 y 2 - h\yi)y 2 . 

(3d) 


Instead of the logistic mode of resource supply, as is the case in MacArthur’s 
model, our model is based on so-called “equable” mode of resource exploita¬ 


tion 


24|, by which the quantities of available resources are held constant by a 


continuous-flow system. According to (l3al) and (|3bl) . a constant concentration 
of jth resource (j = 1,2) flows into a defined volume with the rate pj while 
unused resource flows out with the specific rate qj , in much the same manner 


as in a chemostat 


25| . In natural conditions, the equable modes of feeding, for 


instance, can be found on the first trophic level of ecosystem, among autotrophs. 

We assume that the consumers’ functional response is linear. Endowing 
the consumers with nonlinear (and even distinct from one another) functional 
responses seems a premature complication of the model. Such an extension of 
MacArthur’s generic model, as shown by Abrams and Holt j^dj], may lead to 
several notable modes of behavior, including coexistence via periodic cycling. 
However, the resulting coupled oscillations in their model are rather entrained 
than “awakened”, because one of the two involved species, when unlinked, is 
able to oscillate. 
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Besides, in (l3cl) and (l3dl) intraspecific competition strength di (i = 1,2) 
measures direct interference of individuals within ith consumer population with 
each other resulting in an additional per capita loss rate diyi ; interspecific com¬ 
petition strength h s (s = 1,2; s ^ i) quantifies direct interference effect from 
sth consumer on ith consumer resulting in an additional per capita loss rate, 
h s y s , of the latter. 

Equations ([3]) contain two important assumptions. First, they assume that 
the resources are noninteractive. On higher trophic levels, however, resources 
may interact and the possibility of competition among the resources was origi¬ 
nally pointed out by Levine[27] and empirically confirmed by Lynch [28]. Since 


then, a whole series of theoretical papers (based on MacArthur’s equations) have 
been published on two-predator, two-prey systems with interference competition 
between two self-reproducting prey species 


29 
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As seen from (15al) and there is no intraspecific interference competition 
within the resource populations, in distinction to MacArthur’s model. Yet the 
resource abundance would remain finite even in the absence of the consumer. 

The second assumption of our equations is that the consumers interact only 
directly, through interference competition. They cannot compete trophically, 
through their use of resources, as each consumer specializes on one resource 
only. In the models of pure trophic competition^ like MacArthur’s logistic- 


supply model 


31] and equable-supply model 24 32|, each predator is allowed 


to feed on both prey. 

Intraspecific interference competition is allowed within the consumers as 
well. Owing to this assumption, a consumer would remain bounded even though 
the abundance of the associated resource happened to be constant. 

The novelty of model ([3]) is that it considers time hierarchy of MacArthur’s 
consumer-resource equations to be reversed by assuming dynamics of the con¬ 
sumers to be much faster than that of the involved resources and articulates 
the importance of direct competition mechanisms within the framework of this 
assumption. 

The distinguishing features of aforementioned versions of generic MacArthur’s 
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Figure 1: Variations on a theme of MacArtur’s model as proposed by different authors. C\ 
and C 2 are consumers, R\ and R 2 are resources. Solid arrows indicate predation, dashed 


mean interference, (a) 


0 , 0 ; (b) 00 ; (c) 00 ; (d) 0 ]; 


(e) present paper. 


model, including the system under consideration, are summarized diagrammat- 
ically on Fig. [Tj 

Upon the scaling u\ = CiW\X\/bi — 1, u 2 = c 2 W 2 X 2 /b 2 — 1, v\ = Ciyi/qi, 
V2 = 022/2/92, 7 i = CipiWi/hqi - 1, 72 = C 2 P 2 W 2 /b 2 q 2 - 1, <$i = di<?i/&ici, 
S 2 = d 2 q 2 /b 2 c 2 , = hiqi/b 2 ci, x 2 = h 2 q 2 /bic 2 , /3 = qi/q 2 , £i = qi/h, 

£2 = q 2 /b 2 , and t' = q\t, equations (J3]) take the following nondimensional form; 

Ui = 71 - UiUi - Ui - Vi, 
pu 2 = 72 — U 2 V 2 — U2 — V 2, 

(4) 

£\Vi = (wi — (5iUi — X 2 V 2 )Vl 1 

e 2 V 2 = ( u 2 — S 2 v 2 — XlVi)v 2 - 

Note that in (]4]) dots mean differentiation with respect to nondimensional “slow” 
timescale variable t' measured in units of the resource lifetime 1 jq \, as defined 
by the chosen scaling. 

The parameters /3 _1 , and £^ 1 reflect the rapidity of the dynamics of u 2 , 
Vi and v 2 with reference to that of u\. It is assumed that /3 = 0(1), £i ,£2 1 

and 81,82 < 1 . 

For the sake of simplicity but without any loss of generality, we set /3 = 1, 
£1 = £ 2 = £ and Si = 82 = 8 , and also drop the prime at t. At the same time, 
we retain resource income rates, 71 and 72 , and coupling strengths, n\ and x 2 , 
as free control parameters of the model. Eventually model equations take the 
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form 


«i = 7i - (ui + l)«i - ui, 

U 2 = 72 — (U2 + 1)«2 — U2, 

( 5 ) 

SVi = (ui — <5x71 — X 2 V 2 )l’i, 

£V 2 = ( u 2 — 5 v 2 — XlVi)v 2 . 

It should be mentioned that being proportional to its dimensional prototype, 
Vi directly represents population density of consumer species and is always non¬ 
negative. Quantity u t , however, is not an abundance of resource in the true 
sense of the word. It is rather an affine transformation of x. t done for reasons 
of mathematical convenience. Unlike a purely linear transformation, an affine 
map does not preserve the zero point, so in © Ui = — 1 corresponds to zero 
level of ith resource in reality. Nevertheless, from here on we will apply the 
term “resource” to Ui for brevity. 


3. Analysis and implications 


When x 2 = 0, the consumers are independent. An uncoupled consumer- 
resource system obeys the equations 


u = 7 — (it + l)v — u, 
£V = (u — Sv)v, 

which have a unique stable positive steady state: 


( 6 ) 


u= UV 1 + (47 + 2 + <5)<5- 1-5] = 7 d + 0(S 2 ), 

,_ ( 7 ) 

v = [y/l + (4y + 2 + 5)5 — 1 — <5] = 7 — 7(7 + 1)<5 + 0(S 2 ). 

Equilibrium 0 is a node (intrinsically nonoscillatory steady state) for e = o(S 2 ). 
In the subsequent discussion we assume that this condition is fulfilled. 

Physically feasible equilibria, (ui,u 2 ,vi,v 2 ), of (0 are those for which v \, v 2 ^ 
0. We denote the interior fixed point by F\ 2 = {ui,u 2 ,Vi,v 2 ), where the sub¬ 
scripts at “F” stand for the consumers. Lack of a certain index at a boundary 
fixed point means that the consumer concerned is not present (extinct). Thus 
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Table 1: Existence and stability conditions of nonnegative equilibria in system ©. 


Equilibrium 

Existence 

Stability 

F 

Always 

Never 

Fx 

Always 

72/71 < *7 

f 2 

Always 

72/71 > 1/^2 

F 12 

l/x 2 < 72/71 < *7 

for yt \, >C 2 > 1 (strong coupling) 

< 72/71 < 1/^2 

for x\,*C 2 < 1 (weak coupling) 

Never 

n-i = o(e 1 / 2 ) 


Fi = (ui, M 2 , wi, 0) and F 2 = (u\, U 2 , 0 , 1 * 2 ) designate either of one-consumer 
equilibria corresponding to dominance, while F = (ui,U2,0, 0) means both con¬ 
sumers having been washed out. 

Model (J3J) has four feasible steady states. To 0(1) for small S 


F : 

Ui 

= 71J u 2 = 72, 

Vi = 0, 

V2 = 0 

El : 

Ul 

= 0, U 2 = 72, 

vi = 71, 

O 

II 

, CN 

1^ 

F 2 : 

Ui 

= 71> u 2 = 0, 

Vi = 0, 

V2 = 72: 

F\2 '■ 

Ul 

*77i - ^272 - 

+ 1 ± R 

2(x 1 

-1) 

5 


— 

-Mill + ^272 

— X'l }<2 + 1 ± R 


2(x 2 - 1) 
Vl=U2/*l, V 2 = u 1 /x 2 , 
where R = [(>fi 7 i — ^272 — + l) 2 

+ 4 x 2 (>q - l)(xi7i - 72)] 1/2 . 


( 8 ) 


Existence and stability conditions of equilibria © are summarized in Ta¬ 
ble [I] The model reveals qualitatively different behavior at strong and weak 
coupling between consumer-resource pairs. For lack of space we are not able 
to discuss the case of strong coupling in detail and restrict ourselves to a brief 
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comment. If coupling is strong, any static coexistence of competing consumers 
is not possible: one of the consumers wins and completely dominates. Intense 
competition makes possible bistability of boundary equilibria, as evident from 
Table [1] When both F\ and F 2 are stable with an unstable coexistence steady 
state F 12 , the system being studied is able to exhibit a hysteresis effect. In the 
case of bistability the winner is determined by the initial conditions. 

Application of Routh-Hurwitz stability criterion yields that Fi 2 is stable for 
fairly small coupling strengths, = o(e 1 / 2 ). However, e is so small, that 

from the practical viewpoint, F\ 2 turns out to be unstable for any moderately 
weak (physically reasonable) coupling. As seen from Table [T] the very existence 
of the interior equilibrium F\ 2 in the case of weak coupling, x 2 < 1, implies 
instability of both boundary fixed points, F± and F 2 . System m happens 
to possess four nonnegative steady states, none of them being stable. The 
instability of F± 2 is through growing oscillations. In such a case, the model 
would thus be expected to have a limit cycle in its four-dimensional phase space 
corresponding to sustained oscillations. For sufficiently weak coupling strengths 
of order (^(e 1 / 2 ), i. e. not too far away from the Hopf bifurcation, this limit 
cycle is small and represents a low-amplitude quasi-harmonic periodic solution. 
As a practical matter, the range of such an infinitesimally weak coupling is of 
less concern to us than is the range of far more feasible not-too-weak coupling, 
corresponding to well-developed substantially nonlinear oscillations. 

By the assumption, 0 < e <C 1, meaning that system © is singularly per¬ 
turbed. The slow variables are resources, u\ and u 2 , and the fast variables 
are consumers, v\ and 1)2 ■ The standard practice of reducing such systems is 
multiple-scale analysis 33] whereby fast variables are adiabatically eliminated. 


One has to establish the validity of the adiabatic elimination in each specific 
case. In particular, Tikhonov’s theorem requires quasi-steady states of the fast 
equations to be stable. 

To decompose the full system (J5J into fast and slow subsystems, introduce 
fast time variable r = t/e. Now rescale © by replacing t with re and, after 
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taking e = 0 , it becomes 


u'l =u' 2 =0, 

v[ = (mi - dvi - x 2 v 2 )vi, (9) 

v 2 = ( m 2 — Sv 2 — XlVi)v2, 

where prime means differentiation with respect to r. This is the fast subsystem, 
where u,\ and 112 are replaced by their initial values and treated as parameters. 
It yields the inner solution, valid for t = 0(e). 

Setting e = 0 in © leads to the slow subsystem 


ill = 7i — (ui + 1)^1 - mi, 

(IPa) 

U2 = 72 — (U2 + 1)V2 ~ M 2 , 

(10b) 

P = (Ui - Svi — X 2 V2)Vi, 

(10c) 

0 = (U2 — — XlVi)V2, 

(10d) 


which produces the outer solution, valid for t = 0(1). In this singular limit as 
£ —y 0, the subsystem defines a slow flow on the surface (slow manifold) given 
by (11 Pel) and (llOdl) . Outer solution is valid for those u\ and u 2 , for which the 
quasi-steady states of the fast subsystem © are stable. 

We anticipate the dynamics of the full system ([5]) in its four-dimensional 
phase space (iti, it 2 , iq, u 2 ) to consist of two typical motions: quickly approaching 
the slow manifold (11 Pel) and (llPd|l . and slowly sliding over it until a leave point 
(where the solution disappears) is reached. After that, the representing point 
may possibly jump to another local solution of (II Pci) and (IlPdl) . 

Thus, we ought to find all quasi-steady states of the fast subsystem ([9]), 
map the domains of their stability onto the slow phase plane (zti,u 2 ), and then 
investigate the dynamics of the slow subsystem m with piecewise continuous 
functions. 

The fast subsystem ([9]), which is nothing but the classical LVG model, has 
four quasi-steady states—three boundary and one interior—denoted by Q (the 
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0 



0 ri (5 72*2 7i 

Ul 


Figure 2: The limit cycle of system © projected onto phase plane of the slow variables 
(resources u\ and 112 ). The direction of motion is indicated by an arrow. Consumer 1 is in a 
stable nonzero quasi-equilibrium for any combination of resources below the line h\u\ — 5u2 = 
0, consumer 2—above the line 5u\ — K 2 U 2 = 0. For better appearance we perform scaling 
Ui —> arsinh(wi/7j<5) (i = 1,2). 


slow variables are deemed to be frozen): 


Q : v i=0, V 2 = 0; 

Q 1 : vi= ui/5, v 2 = 0; 


Q 2 '■ 
Ql2 ■ 


vi =0, v 2 = u 2 /5 ; 

X 2 U 2 — Sui _ X\Ui — 8U2 
Vl r? ' ^2 r~2 

X1X2 — 0 X1X2 — 0 


(lla) 

(lib) 

( 11 c) 

(lid) 


Quasi-equilibria Q\ and Q 2 are stable nodes respectively for 6112 < *c\U\ and 
Sui < X 2 U 2 ', Q and Q 12 are always unstable. 

Consider the plane of resources shown on Fig. [2] Let initially Qi be sta¬ 
ble and Q 2 unstable with consumer 1 dominating. This corresponds to slow 
variables u\ and 112 being somewhere below the line Su\ — X 2 U 2 = 0. The dy¬ 
namics of the resources (treated as bifurcation parameters in reference to the 
consumers) is described by a system of two independent equations 


ui = 7i - [('»i + !)/£ + l]ui, 
U 2 = l 2 ~- U 2 , 


( 12 ) 


12 











which, in view of ([ITbJl . is a oiecewise version of the slow subsystem m- System 
m has a stable steady state 


^ — ^[\/l + (471 + 2 + S)S — 1 — 5\ — 7i(5 + 0(S 2 ), 


-(1) 

= 72 . 


(13) 


It is marked by an open circle in the upper left corner of Fig. [2] While heading 
to m, the trajectory crosses the line Sui — X 2 U 2 = 0 and enters the domain of 
bistability of both Q\ and Q 2 - However the dominance of consumer 1 persists. 

Note that in (Mb the variable u\ is faster than 112 due to small 6 . Clearly, 
the representing point must have relaxed to the vertical line u\ = ~ 7i<5 

well before approaching the horizontal line u 2 = = 72 - 

Eventually the trajectory has to cross the line xqui — Su 2 = 0. As soon 
as this has happened, node Q\ of the fast subsystem © will be absorbed by 
saddle Qi 2 - Consumer 1 rapidly washes out, and Q 2 becomes the only stable 
quasi-equilibrium, with consumer 2 dominating. 

In terms of the four-dimensional phase space of full system ([5]), the repre¬ 
senting point is now in the other stable branch of the slow manifold given by 
(ll()c|) and (llOdl) . The motion over this alternative branch obeys the piecewise 
subsystem 


Ul = 7l - Ml, 


(14) 


U2 = 72 — [(«2 + l)/5 + 1]m 2 

with the initial conditions Mi(0) = 71 A and ^ 2 ( 0 ) = 7 i>q. The dynamics of (1141) 
is basically similar to that of (fl2l) analyzed above. The variable 112 compara¬ 
tively rapidly relaxes to = 72 A-F C>(<5 2 ); the variable u\ slowly grows toward 
u\ = 71 . When u\ has crossed the level 72 ^ 2 , node Q 2 would be absorbed by 
saddle Qi 2 - The system returns to the first branch of the slow manifold, and 
thereby the oscillatory cycle gets closed. 


4. Results and discussion 

As seen from Fig.[2l the two coupled consumer-resource communities execute 
self-sustained synchronous antiphase-locked oscillations. Fig.[3]shows the results 
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Figure 3: Time profiles of self-sustained relaxation oscillations experienced by coupled 
consumer-resource pairs, (a) resource 1 , (b) consumer 1 , (c) resource 2 , and (d) consumer 
2. The numerical values of the parameters are e — 0.73x10“ 3 , 71 = 1.2, 72 = 1, xi = 0.5, 
X 2 = 0.8, and S = 0.1. 


of numerical integration of ©. The resources demonstrate sawtooth periodic 
pulses. The oscillation range for the resource levels remains finite and, what 
is important, it does not depend on the intraspecific interference parameter 5. 
The consumers change periodically between extinctioij^ and respective constant 
levels 71 and 72 . Very brief transient from zero to flat nonzero level within each 
cycle is accompanied by a highly pronounced spiky overshoot. The magnitude 
of the spike tends to infinity as S —> 0, in view of (lllbl) and (lllcl) . 

One may distinguish four parts within the period of synchronous oscillations: 

1) Consumer 1 is essentially zero, while consumer 2 is approximately equal to 
its uncoupled steady-state value, 72 . Resource 1 increases due to its constant 
inflow until it overcomes losses for consumer 1 ; 

2) With a sufficient resource stock, consumer 1 now emerges. The population 


2 Actually, the exact solution to full system m yields nonzero v\ and V 2 at any time, even 
though they may take on very small values, so that ln(u) oc — e~ 1 . It should be emphasized 
that the competing consumers periodically “die out” only within the framework of a multiple- 
scale technique’s approximation being used. This is merely convenient idealization. 
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1 exhibits a spike due to the fast time scale of the consumer equations. The 
sharp increase in population saturates the available resource level, so resource 
1 drops. Cross-losses cause consumer 2 to wash out; 

3) Quantities V\ and u\ relax to their equilibrium values, as if there were only 
one uncoupled consumer-resource pair. Consumer 2 is essentially zero. Re¬ 
source 2 is increasing, like resource 1 did in part 1; 

4) Resource 2 surpasses the losses, consumer 2 emerges and the subsequent 
cross-losses cause consumer 1 to wash out. The spiking consumer 2 also 
causes a substantial decrease in the available stock of the associated resource. 
The sequence begins again. 

The essential feature of the model is that when one consumer is very scarce, the 
whole coupled system behaves like an uncoupled consumer-resource pair ([6]). 

Presenting his famous model Smale remarked that “it is more difficult to 
reduce the number of chemicals to two or even three” jf], p. 26]. As distinct 
from Smale’s example, the bilinear coupling in our case makes self-sustained 
synchronous oscillations possible for just two variables. 

As we have seen, phase trajectory of the system constantly moves from the 
neighborhood of unstable boundary equilibrium F\ where only consumer 1 is 
present, to the neighborhood of F -2 where consumer 2 completely dominates, 
back to Fi , and so on in cyclic alternation. This kind of trajectory was termed 
“heteroclinic cycle” by Kir linger 21,. A heteroclinic. cycle occurs when the 


outflow (unstable manifold) from one saddle point is directly connected to the 
inflow (stable manifold) of another saddle point, and vice versa. It is closely 
related to another notion of the nonlinear dynamics, a homoclinic cycle, which 
emerges when the unstable and the stable manifolds of the same saddle coincide 
and form a closed loop. 

Homo- and heteroclinic cycles are not robust structures in the sense that 
infinitesimally small change of system parameters destroy them. However in 
the practical sense, any limit cycle passing in close proximity to saddle points 
will be indistinguishable from a heteroclinic cycle (Fig. [4]). The only difference is 
strict periodicity, although the period of the limit cycle in a neighborhood of the 
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Figure 4: A 3D-projection of the limit cycle in system © for parameters chosen in a neigh¬ 
borhood of the heteroclinic cycle. 


heteroclinic cycle may be long. Besides, at the threshold of homo-/heteroclinic 
bifurcation the period is susceptible to external noise. 

In the context of our model, as coupling becomes stronger, the stable limit 
cycle swells and passes closer and closer to boundary fixed points which are 
node-saddles. Depending on the interplay between the parameters, eventually 
it may bang into one or both of these equilibria creating either a homoclinic or 
heteroclinic cycle, respectively. This corresponds to 72/71 = x'i and 71/72 = >^2- 
On further increasing the coupling, the saddle connection breaks and the loop 
is destroyed. 

It is worth noting that heteroclinic cycles were first found by May and 
Leonard £ in a classical LVG system with competing three species. How¬ 
ever in their model the cycle is not truly periodic: as time goes on, the system 
tends to stay in the neighborhood of any one boundary equilibrium ever longer, 
so that the “total time spent in completing one cycle is likewise proportional to 
the length of time the system has been running.” Moreover, May and Leonard 
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state that “the phenomenon clearly requires at least three competitors, which 
is why it cannot occur in models with two competitors.” This statement is 
echoed by Vandermeer [35] who extended their theory on higher dimensions: 


“It appears to be the case that all cases of an odd number of species follow this 
basic pattern, whereas all cases of even number of species result in extinction 
of half of the components, leaving the other half living independently at their 
carrying capacities.” In view of our results, the above conclusion is by far and 
away true providing one stays within the framework of classical LVG equations, 
which in fact implies a high rapidity of the resource dynamics. In our model 
of just two competitors the slowness of the resource relative to the consumer is 
essential for the oscillations to occur, because it provides the necessary inertia 
to the system. 

The feasibility of our model is tightly bound to justification of the adopted 
time hierarchy in system (|5|). In ecosystems, the most common case is rapid 
consumption of food by species. However it seems reasonable to propose that 
the model may describe the first level of an ecosystem, at which the consumers 
are autotrophs and the resources are mineral nutrients. The ability to exploit 
different substrates leads to a possibility of stable coexistence of different or¬ 
ganisms descending from a common ancestor. Divergent evolution is just the 
emergence of new species: due to mutations two populations come into be¬ 
ing, sharing the same genetic code but having proteins able to process different 
substrates. Providing the environmental conditions are quite stable on the evo¬ 
lutionary timescale, the inflows of inorganic substrates from the surroundings 
may be considered constant and the washout time of a substrate may occur 
much longer than the life expectancy of a species (recall the definition e = q/b). 

From a non-ecological perspective, by and large similar relationships can 
be found in coupled longitudinal modes of laser with second harmonic genera¬ 


tion. Baer 123 experimentally observed antiphase oscillations of two (and more) 
modes in a multimode neodymium-doped yttrium aluminum garnet (Nd 3 + :Y 3 Ali 50 i 2 ) 
laser with an intracavity potassium titanyl phosphate (KTiOPCH) frequency¬ 
doubling crystal. He also proposed a model for that phenomenon in terms of 
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the coupled rate equations. In original designations, the equations for the pop¬ 
ulation inversions, G\ and G 2 , and the intensities, I\ and / 2 , are as follows: 


T f Gi=G°i-(PiIi+/3i 2 l2 + l)Gi, 

(15a) 

T f G 2 =G° 2 - (p 2 h + fah + 1 )G 2 , 

(15b) 

t c Ii = (G 1 — ai — eli — 2eI 2 )Ii , 

(15c) 

T c i2 = (G2 — 0-2 — el 2 — 2eli)l2, 

(15d) 


where Tf and r c are the fluorescence time and cavity round trip time, respec¬ 
tively; Pi and @2 are self-saturation parameters, which determine how strongly 
the corresponding intensity depletes the available gain; /3 i 2 and fa are cross¬ 
saturation parameters; Gf and G 2 are the small signal gains (pump parameters); 
oc\ and a .2 are the cavity losses for the respective modes; and e is the nonlinear 
coupling coefficient due to the presence of the intracavity doubling crystal. In 
(115cl) and (115dl) . the second-order terms if, if and I1I2 account for the loss in 
intensity of the fundamental frequencies through second harmonic generation 
and sum-frequency generation. 

It is apparent at a glance that m represents coupled consumer-resource 
equations, where the intensities play the role of consumers, while the population 
inversions act as resources. Structurally, Baer’s system is very nearly the same 
as system (O being discussed. Even the time hierarchy is similar: r c /ry = 
(0.5 ns)/(0.24 ms) = 0(1O~ 6 ), i. e. the “resources” change much slower than the 
“consumers”. There are three distinctions, however. First, in (115cl) and (II Sell) , 
intra- and interspecific interference parameters are not independent, both of 
them being proportional to the coupling strength. Second, in (I15al) and (I15bl) . 
the two modes are allowed to “compete” for the active medium, so the consumer- 
resource pairs turn out to be further linked trophically. Third, when uncoupled, 
the steady state of each mode is stable focus, not stable node as is the case in 
our model. In other words, uncoupled consumer-resource pairs in Baer’s system 
are weakly damped oscillators with the intrinsic period proportional to ^ /T c Tf. 
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Baer performed numerical integration of (1151) for different coupling strength 
£ and with many different initial conditions. He revealed that as e decreases, 
the mode-coupling oscillation period decreases tending to the period of intrinsic 
oscillations. If e is decreased further, the oscillations cease, and the system 
becomes stable. With a large e, the oscillation period becomes quite long, and 
each mode appears to reach a stable intensity value before abruptly switching 
off (cf. the results of our analysis!). The numerical solutions correctly predicted 
that the two modes tend to pulse on and off out of phase with each other. 

Subsequently, Erneaux and Glorieux js7, pp. 318-325] reduced (fl5l) to the 
equations for coupled quasi-conservative oscillators and proved the existence 
of stable antiphase periodic solution in the case of the modes with identical 
parameters G°, a and /?. However that result has to do with the onset of low- 
amplitude quasi-harmonic oscillations. Unlike their study, our approach deals 
with well-developed high-magnitude essentially nonlinear oscillations. 

Interesting issues concern how the outcome of exploitation and interference 
is altered when the mode of resource supply is not constant, or when the in¬ 
terspecific interference is not necessarily mutually costly (i. e. each consumer 
suffers a net reduction in per capita growth rate via interference from, but can 
gain an increase in growth rate via interference on, the other consumer) 38]. 


Investigations of these possibilities may constitute a future direction for work 
on the model. 


5. Conclusions 

We considered a model of two consumer-resource pairs linked by interspecific 
interference competition. When uncoupled, an individual consumer-resource 
pair has a unique stable steady state and does not admit periodic solutions. If 
intraspecific interference within the species is strong enough, the equilibrium is 
nonoscillatory. 

When coupling is moderately weak, the model reveals low-frequency an¬ 
tiphase relaxation oscillations. The consumers cannot coexist even dynami- 
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cally: in each of two periodically alternating states one consumer completely 
dominates and the other is on the verge of extinction. The most intriguing 
feature of the model is that each of the involved consumer-resource pairs taken 
separately does not oscillate; both communities are completely quiescent, how¬ 
ever, in interaction, when coupled in a nonlinear way, the resulting system turns 
into a relaxation oscillator. 
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